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In order to improve the control accuracy of the robot manipulator, the sliding 
mode control combined with the adaptive neural network (ANNSMC) is 
proposed. Sliding mode control (SMC) is a nonlinear control recognized for 
its efficiency, easy tuning and implementation, accuracy and robustness. 
However, higher amplitude of chattering is produced due to the higher 
switching gain to handle the large uncertainties. For the purpose of reducing 
this gain, the uncertain parts of the system are estimated using neural 
network (NN) with on-line training using back propagation (BP) technique. 
The results of the online interconnection weights between the input and the 
hidden layers and between the hidden and the output layers are injected 
offline in order to improve the network performance in term of the 
convergence speed. In order to reduce the response time caused by the 
online training, the obtained output and input weights are updated using the 
adaptive laws derived from the Lyapunov stability approach the proposed 
control ANNSMC has improved the convergence speed with 41.13% for the 
first link and 40.15% for the second link comparing to NNSMC. The 
simulation result illustrates the performance of the proposed approach by 
using MATLAB and the control action suggested did not manifest any 
chattering behavior. 
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1. INTRODUCTION 


The motion control design for two links robot is an intense area of research that has fascinate 
considerable focus, it’s a big challenge because of the nonlinearities of the system and the uncertainties of 
their parameters. Furthermore, the right dynamic models are nearly hard to obtain because the system is 
described by a nominal model with major uncertainties, to name a few: external disturbance, internal friction, 
and payload parameter. Several strategies have been proposed to cope with parameter uncertainties including 
neural network (NN) based controls [1]-[22], neural adaptive proportional-integral-derivative (PID) control 
[11], the sliding mode control (SMC) [23]—[31], SMC and PID controller tuned by whale optimizer algorithm 
(WOA) [29], nonlinear model predictive control tuned by NN [32]-[34] , the adaptive fuzzy control [35]— 
[37], the particle swarm optimization (PSO) combined with NN and SMC [38], the combined PID and SMC 
[30] were, the self-tuning fuzzy PID-nonsingular fast terminal sliding mode control for robust fault tolerant 
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control of robot manipulators is studied. Sliding mode control (SMC) is an important robust control, it 
provides a systematic approach to the problem of uncertainties, bounded external disturbance and 
nonlinearities [5], [30], [32]-[37]. 

Unfortunately, there is an undesirable chattering in the control effort and it involves extremely high 
control activity. It is worth mentioning that eliminating this chattering problem several solution are proposed 
to perform properly as the boundary layer solution [23]-[25]. However, this method in only efficient for small 
uncertainties. In case of large uncertainties, the unknown parts of the two-links robot are estimated by using a 
neural network architecture, this approach leads to fewer uncertainties and a lower switching gain as result. 

The back propagation technique [4], [38], [39] is used to train the weights of the neural network on- 
line. The most important task is to ensure the system’s asymptotical stability and the tracking error’s 
convergence to zero. In this study, the online input and output weights are injected and offline changed using 
the adaptive rules obtained from the Lyapunov stability theorem [40], therefore the asymptotical stability of 
the system is ensured, and also the speed of convergence is widely improved and illustrated in simulation 
results. In addition, we used integral squared error ISE [41] to validate numerically the theoretical statements, 
and that gives outstanding results. 

This paper is organized as follows. The dynamic model of two links robot systems is introduced in 
section 2. Section 3 covers the neural network design technique; outlines sliding mode control and illustrates 
the adaptive laws for modifying the input and output weights obtained from the Lyapunov synthesis approach 
and the stability proof is presented. Section 4 displays the simulation results obtained by using a two-link 
robot manipulator. Finally, the conclusion is addressed in section 5. 


2. SYSTEM MODEL OF THE TWO LINKS ROBOT 
Consider the dynamic model of the two links robot written as (1) [7]. 


M(Q)4 + Vn(a Ga + Fq + f-(q) + G(q) + tq =T () 


With g,q,q € R’represent respectively the joints position, velocity and acceleration. M(q) € R**? is a 
symmetric, positive definite inertia matrix, V,,(q,q) € °*? is the centripetal and coriolis matrix, F € 92? 


denotes the viscous friction coefficients, f,(q) € R? is the coefficients of coulomb friction, G(q) € R? is the 


vector of gravitation, Tg € $2” represents unmodeled dynamics and bounded unknown disturbances, t € SR” 
is the input torque vector applied to the servo motor. 


V, V, M,,M 
V qQ)j= m1 na M =. 11 a 
m4 4) Vinzi_ Vm22 @) M21 M2 
m1? 4mz1? 5 : m1? m1? c m1? 
My, = 3 +—J— + ml cos(qz) 3 Miz = M2, = 3 + 5 cos(qz); Mz2 = 5 
With q = [%1 2]" the positions; 1 = 1, = l,: the lengths; m1, m,: the masses. 
ml? . , ml? .. : . 
Vinit = 42 sin(qz2), Viz = — 7 (G2 + 41) sin( gz) 
mygl mzgl 
ml2. ., m2 gl cos(q,) + 2 cos(q,) + 2 cos(q, + G2) 
Vinzi = arr by sin(q2) Vm2z2 = 9;G(q) = mgl 
— c0s(h + q2) 


With g is the gravitation term. 


fards + ky sgn( qi) 


H=Fqt+f-(q) = faa + ky sgn Gz) 


System model described by its state space is as (2). 


xX, = Xz 
Xq = fi (X) + G11 1, X%3)Uy + G12 (%1,%3)U2 + G1 (x, t) (2) 
X3 = Xq4 


X4 = fo(X) + Gor (%1, X3)Uy + G2 (%1,%3) U2 + % (x, t) 
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Where: ¢,(x,t)and (x,t) the unknown parts, tT=U=[U1 U2]? =[% T2]’q=[% %3]",g= 
911%, X3) seo ath 


G2i%1,%X3)  J22(%1,X3) 
911 (%1,X3) > Nand gz2(x1,x3) > 0 and f@) =[fi@®) foO]" =—-M7*{-V,,[%2 *4]" — H-G}. 


a a 
[X2 X4).G@ =[%2 x4), x= [eaxxgxq]" and M7! = g(x,,x3) = ( 


3. ADAPTIVE NEURAL NETWORK SLIDING MODE CONTROL DESIGN 
3.1. Controller design 

This study examines the trajectory tracking problem of the robot manipulator mentioned previously. 
The challenge of trajectory tracking control is to construct a robot controller based on any of x(0)eR”, 
x(t) — xg(t) > Oas t > oo. Where x(t) and xq(t) are respectively the trajectory and the desired trajectory of 
a robot. In the following, some variables are defined as (3). 


§@,) P= [ht PPP &@0] (3) 
The system's output tracking error can be described as (4). 
€e=q-qq =[%17— %1a X3 — X3q]7 


€=q-qa= [% — %1a x3 ~ da)" 
€=G-Ga=[%i-—hiaq ¥%3 — Xa)" (4) 


With x, and x4,q are the desired output. r = 2 is the relative degree of the system (2), and the sliding surface 
is characterized as (5). 


S=é+ fe (>) 


Bu =O 
0 Boz 


B is designed diagonal matrix as B = ( ) The selection of # must satisfy the following Hurwitz 


polynomial as (6). 


s@ + Bis =0 
s + Boas = 0 (6) 


With s = <o The sliding variable derivative is (7). 
$= fy + §(t) + gu - (524) + Be a) 
4d 
The following is the robust adaptive controller that was used (8). 


U= Un + Us (8) 


with u, = —ksign(S). The control law that respect (8) is (9), 


w= 951(x) (~falx) + (G24) — Be) — ksign(s) 0) 
where sign( )is the sign function given by: 
lifS > 0 
sign(S) = j0ifS = 0 
—1ifS < 0. 


In order to compensate the uncertainties, the positive switching gain k is designed as B < k. With B the 
upper bound of uncertainties given by lz (x, t)|| < B. The boundary layer method may be utilized to remove 
the chattering effect induced by the discontinuous control law. The control is now as (10), 
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x ‘ 
w= 951(x) (~falx) + (G24) — Be) — ksat(s) (10) 
with the saturation function is described by (11), 
_ (S/dif[|S|| < 6 
oe esas sense (it) 


with 6 is the boundary layer thickness. 

This approach is suitable for systems with minor uncertainty. We suggest using NN to estimate the 
uncertain terms of the system presented in (1) for large uncertain systems, in order to keep system 
uncertainties manageable. Let denote the prediction of unknown nonlinear functions parts as: & (x,t) = 


[2@t) €(x,0)]’, with as (12): 
e(x,t) = F(x, t) — §@,t) (12) 


and: lz (x, t)|| < e* where é* is the network prediction errors’ upper bound. 
- Theorem 1: Consider the robot manipulator stated in (1) in the context of significant uncertainty. If the 
system control is programmed as (13): 


u= gi'(a) (— (C2) + £8) + (228) — Be ~ ksat(s)) a3) 


with e* <k and é (x, t) is predicted by the proposed off-line NN structure, then the trajectory tracking 
errors converge to zero in limited time. 
- Proof. Consider the potential Lyapunov function: V = ) Ts then V = S’S replacing the expression of S$ 


given in (7) we have: 
V = ST + Et) + out (FA) +78 + 88) 
4d 


by substituting the expression of 4 provided in the theorem we get: 


V = S™(E(x,t) — €(x, t) — ksat(S)) = STe(x, t) — kSTsat(S) < ||S™||lle(x, t)|| — kSTsat(s) 
< ||STlle* — kS7sat(S) 


by selectinge’ < k, with k s a minor gain that is just responsible for compensating network faults forecast, 
we have: if ||S|| = 6 then we take sat(S) = sign(S) for anyd > 0 and the function V = (e* — k)|IS|| < 


Ss : tic La 
0. However, we take sat(S) = 52s continuous function in a boundary layer (a small 6 -vicinity of the 


origin), as a result the system trajectories are restricted to a sliding mode manifold boundary layer. S = 0. 
Basing on (10), (8) and (4) we can write as (14). 


é = Aé + BC (x,t) — Bon(w)us (14) 


The neural controller aim is to drive the system output to match the defined intended reference trajectory 
as closely as feasible. The suggested neural network's design technique will be described in the following 
paragraph. 


3.2. Neural network representation 

We investigate a NN with two layers of tunable weights in this study [8]. The architecture used is 
illustrated in Figure 1 with one hidden layer, where X is the state input variables: the joints position, and the 
output variables are the unknown parts given by (2): y, = & ee t) and y, =, (ay t). 

yy = Wi a(W;"x) k = 1,2 Where o()is the activation function for hidden-layer considered as a 


sigmoid o(s) = — function given by Wy = [Wi Wi2Wen]" and W; = [W,1W;2 Win are the connectivity 


e7s 
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weights between the hidden and output layers, as well as between the input and hidden layers. The actual 
output Vax (x) which is the difference between the actual and nominal functions is (15). 


Vak (x) = yx (x) + e(x) (15) 


Where the approximation error of NN is e(x). The weights of the network are modified during the offline 
implementation. The technique used is based on the BP algorithm, it’s a gradient-descent method and widely 
used. The back propagation technique is used to determine the necessary modifications after randomly 
selecting the network weights. The algorithm may be broken down into five steps: 

i) Forward propagation of operating signal 

i1) Calculate loss function 

iii) Back propagation of error signal 

iv) Calculate gradients 

v) Weight updates 

When the value of the error function becomes sufficiently minimal, the algorithm is terminated. It worth 
mentioning that BP algorithm converges slowly. In order to improve the speed of convergence, the adaptive 
principles obtained from the Lyapunov stability theorem is used to update online the obtained inputs and 
outputs weights. The detailed steps of the proposed algorithm are illustrated in [40]. 


x; ry = €(x,7) 
*'3 
Xy Le etek F100) 


Input Layer Output Layer 


Hidden Layer 


Figure 1. A multilayer neural network's design for predicting unknown components 


3.3. Implementation of adaptation laws 

The network weights are adjusted using the hybrid BP algorithm which takes important time to have 
a result, to deal with this time response weights are adjusted offline. In this case the output of ANN with 5 
hidden nodes can be presented by (16). 


C(x, t) = Wy, * 0% (x, W;) (16) 


The parameters W; and W; need to be adjusted further for the purpose to minimize approximation errors. The 
adaptive rules for them were developed as (17) [23], 


i = —n,0/ B’Pé (17) 


W; = —n2x" BT Pé 


where 7, and 7 are constants that are always positive. P is the positive and symmetric definite matrix that 
corresponds to: 


J =—(ATP + PA) (18) 


where the designer selected J as a asymmetric definite matrix. 
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- Theorem 2: Suppose the nonlinear system described by (1). If the adaptive neural control rule mentioned 
in (13) is used with the parameter adaptation laws (17), as a result, the tracking errors converge to zero 
as t > and all signals in the closed-loop system are limited. 

- Proof: Take into consideration the possible Lyapunov function, which is: 


V =-6"Pé+— WW, 
2 21 
The Lyapunov function's derivative is stated as: 
V ==(é"Pé + 6" Pé) + We Wi 
1 
using in (14) we have: 
V = (Ae + BE(x,t) — Bon (x)us)" PE + 67 P(AC + BE(x,t) — Bon (X)us)) + —We Wi 
1: 
applying (16) and (17) we get: 
V => ((Aé + BE (x, t) — Bgn(X)us)"Pé + €TP(AS + BE(x,t) — Bgn(x)us)) + 
— (=mo¢ BT Pe)" Wy 
uh 
V = <6"(ATP + PA)é += (67 (x, t)B" Pe + 6" PBC (x, t)) —> (ul BT gh(x)Pé + 6" PBg,(x)us) — 
é’PB 0,W,, 
P is symmetric, we get: 


y=? 


~ 2 


éT(ATP + PA)é — ul BT gl (x)Peé 
ee ee eee ; 
Vs 3éTJé —>|lus||B" gn @Plléll s 0 
Hence Vis negative semi definite, the signals éand W, are all bounded. The parameters W, and W; 
described in (17) are adjusted using the projection algorithm as (19): 
—n,o'B' Pe if ||W,||< M,or 


ae Mg and 
e’PBo'W, = 0 


W, = 
_ -TpT pp _ ¢rPBOTW ; ee Mz and 
mo BPe—— Tae We Sterppotw, <0 
=x" BT Peif ||W,l|< Mg.or 
if |Will= Moz and 
W, = e™PBx™W; = 0 


e’PBx'W, W,||= Mg. and 
—1px™BT Pe — j ; {| ilI= Mee 


\|w; ||" e PBx™W, <0 (19) 


The utilization of projection algorithm has a good performance on the tracking trajectory and also in the 
control law illustrated in the next section. 


3.4. Integral squared error (ISE) 

In advance of talking about how to configure a controller, we must first define what makes a 
satisfactory response. It is a great challenge to take up. In reality, there are several metrics that may be used 
to compare the quality of regulated replies. The control measure described in this section is integral squared 
error (ISE) which is used to illustrate the success of the designed control ANNSMC. According to [41], the 
mathematical expression for ISE is (20), 


ISE=f e2dt (20) 


where t is the simulation time, and e is the system's output tracking error. 
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4. SIMULATION RESULT 
In this part, we put the suggested control strategy to the test on a two-link robot represented by the 
model (1). The aim of the control is to keep the system tracking the desired angle trajectory:x,g = 7 


=cos( 0.5t) and x34 = (1/3) sin(0.5t). The criteria are thought to be :m, = m, = 1;1, =1, = 1 the initial 


conditions are: q(0) =[x,(0) x3(0)]’ =0 ,q(0)= [x2 (0) x,(0)]' =0. The uncertainties under 
consideration are vector random noise with a value of unity |tg| <1. The following are the parameters 


100 0 
; ‘aa Mg = 1.35, Mgo = 2.n, = 6,2 = 1.5, A=[1 0:0 1], 


B=[1 0;0 1], P=[-50 0;0 -50], the number of hidden nodes is 5, the gain is k, = kz = 1.5. 

The simulation results show position tracking for links 1 and 2, which are depicted respectively in 
Figure 2 and Figure 3, where the reference signal is represented by the dashed line (red), the results using 
adaptive law for weigh (ANNSMC) is represented by solid line (green), the results without adaptive law for 
weigh (NNSMC) is represented by dashed line (blue) and the results using conventional SM is represented by 
dot line (black). According to Figure 2, the gap between the ANNSMC outcomes position of link 1 and the 
reference value is high at first, but the system remains stable, and the system outputs faster converge to the 
intended trajectory than NNSMC, however the SMC result position don’t converge as the ANNSMC result 
position. Similarly, the ANNSMC position of link 2 closely matches the reference signals and quickly than 
NNSMC. The corresponding control torque signals given in Figure 4 and Figure 5 are smooth even in the 
presence of significant uncertainty, there is no oscillatory behavior. 


connected with the controller design: J = ( 


2.5 r 
: | {\ 
X 2.77 


1.5 | ¥ 1.356 
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Figure 2. Angle response x, 4yysmc Using adaptive law, X;yjwsmc angle response without adaptive law 
X1smc angle response with conventional SMC and desired trajectory x1q 
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Figure 3. Angle response x3 ,nysucUsing adaptive law, X3nvsuc angle response without adaptive law x3>5,y 
angle response using conventional SMC and desired trajectory x3 
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Figure 4. Control u, (input torque of join actuator 1) 
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Figure 5. Control uz (input torque of join actuator 2) 
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In addition, the suggested adaptive neural control method's performance evaluation in term of the 
speed of convergence is investigated in Figure 6 and Figure 7, that represent the position tracking for link 1 
and link 2 where the reference signals is represented by dash red line and the results without adaptive law for 
weight is represented by solid line in blue, by comparing the time response made by the results using this 
approach, Figure 2 and Figure 3, and the results without it, Figure 6 and Figure 7, it’s obvious that the time 
response is widely improved. Furthermore, the error response is also reduced by using the proposed control 
and illustrated in Figure 8 and Figure 9. 


Figure 6. Angle response x, yysyc Without adaptive law, desired trajectory x14 
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Figure 7. Angle response X3yysyc Without adaptive law, desired trajectory x34 


— — —e1ANNSMC 


response error of q1 (rad) 


Figure 8. Error response @;yysmc >» C1annsmc and €;syc of angle q,, respectively without and with adaptive 
law and SMC 


— — —e2ANNSMC 


response error of q2 (rad) 


Time(s) 


Figure 9. Error response €zyysmc> C2annsmc aNd €25myc Of angle gz, respectively without and with adaptive 
law and SMC 


The performance comparison of the ANNSMC control in Table 1 shows that it outperforms the 
neural network sliding mode control (without the adaptive law derived based on Lyapunov theory). The 
control system designed to minimize ISE tend to rapidly reduce huge mistakes, this leads to fast responses 
This is seen in Figures 2, 3, 6, and 7. 
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Table 1. Comparative ISE for ANNSMC and SMC control 
NNSMC ANNSMC SMC 


ey e2 1a e2a C15 25 
ISE 6.1398 80.2373 2.8940 47.241 9.6930 108.6011 


5. CONCLUSION 

The robust reference tracking problem for two-link robot manipulators was solved in this study. The 
developed method is a hybrid of the sliding mode control SMC technology and the adaptive neural network. 
ANN used online to extract weights adjusted firstly with BP algorithm. However, these weights are injected 
to be adjusted secondly by the adaptation rule obtained from the Lyapunov stability theorem. The ability of 
ANN to approximate fatly the uncertainties and external interference were assessed through comparison with 
an traditional SMC methods. Simulations conducted on computers of a two-link robot manipulator confirm 
outstanding results and show that the proposed algorithm can quickly approach the desired trajectory and 
effectively suppresses the chattering. Furthermore, the measure of ISE is used in addition to verify the 
theoretical statements achieved. 
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